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Abstract: This paper argues that automated proofs of identities for non-terminating hy- 
pergeometric series are feasible by a combination of Zeilberger's algorithm and asymptotic 
estimates. For two analogues of Saalschiitz' summation formula in the non-terminating 
case this is illustrated. 
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1. Introduction 

The Gosper algorithm and in particular the subsequent Zeilberger algorithm and related 
WZ method have been extremely successful for an approach by computer algebra to identi- 
ties involving ((/^hypergeometric functions (see the book by Petkovsek, Wilf & Zeilberger 
[7] and the references given there). Most of the current theory and all implementations of 
Zeilberger's algorithm remain restricted to the case of identities for terminating hypergeo- 
metric series, while many of the known identities in the literature involve non-terminating 
hypergeometric identities. The book [7] discusses in Chapter 7 some nonterminating cases 
obtained from terminating cases by the WZ method. A very convincing demonstration 
that Zeilberger's algorithmic is suitable for non-terminating cases is given by Gessel [3, 
p. 547]. He demonstrates Gauss' summation formula for the Gaussian hypergeometric se- 
ries of argument 1 by means of a combination of Zeilberger's algorithm and the asymptotic 
estimate 

(1.1) 

(see for instance Olver [6, Ch.4, §5.1] for a proof of (1.1)). 

I will demonstrate in this paper that the method extends to the non-terminating 
generalization of Saalschiitz' summation formula for a terminating Saalschiitzian 3F2 hy- 
pergeometric series of argument 1. This nonterminating case has the form A + B = C, 
where A is a non-terminating Saalschiitzian 3^2(1), B is a quotient of Gamma functions 
times another 3^2(1) and C is a quotient of Gamma functions. Without knowing the 
formula explicitly, one can start with A, apply Zeilberger's algorithm to the terms of A, 
and finally arrive by some limit transition using (1.1) at the desired formula A + B = C. 

It is important to notice that all steps in the proof, both in the Gauss case and in 
the non-terminating Saalschiitz case, may be automated, including the application of the 
asymptotic formula (1.1). I strongly believe that this idea may lead to an algorithmic 
approach to most identities for non-terminating hypergeometric series in literature. For 
the (/-analogues of such identities an algorithmic approach should be feasible as well. 

For an introduction to Zeilberger's and related algorithms I refer to the book [7]. For 
hypergeometric series see Erdelyi [1] and Gasper & Rahman [2]. Although the emphasis 
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in [2] is on the g-case, it also contains some information on the q = 1 case. Anyhow, many 
formulas for q = 1 can be looked up from the corresponding g-case in [2] , by silently taking 
the formal limit for q j 1, and by using that 

(see [5, Appendix B] and references given there). 

The contents of this paper are as follows. Section 2 discusses Gauss' summation 
formula following Gessel [3, p. 547]. Section 3 discusses the case of a non-terminating 
Saalschiitzian 3 F 2 (1). Section 4 gives three other identities related to the one derived in 
section 3, and it is shown how these identitities follow analytically from each other. Finally, 
in section 5, we consider one of the other identities in section 4 (A + B = C with C a 
Gamma quotient and A and B a Gamma quotient times a non-terminating 3^2(1) with one 
of the upper parameters equal to 1). As is demonstrated, it can be derived by computer 
algebra and asymptotics, similarly as for the case of section 3, but the asymptotics is quite 
tricky (maybe interesting in its own right) and not yet suitable to be automated. 



2. Gauss' summation formula 

Gauss' summation formula (Gauss, 1813; see for instance [1, 2.8(46) with proof in §2.1.3]) 
is the non-terminating analogue of the Chu-Vandermonde summation formula. It reads as 
follows: 



a, b 



;i 



k=0 



{a)k(b) k = T(c)T(c-a-b) 
(c) k k\ T(c- a)T(c-b) 1 



(2.1) 



where we assume that Re (c — a — b) > and c ^ {0, —1, —2, . . .} in order to ensure that 
the hypergeometric series on the left-hand side is well-defined and absolutely convergent. 
Gessel [3, §7] shows how to prove this identity by the WZ method, as I will recapitulate 
now. 

Formula (2.1), with c replaced by c + n (n e {0, 1,2,.. .}), can be written as 



^2f(n,k) = s(n), 



(2.2) 



k=0 



where 



f(n,k) := 



(a)k (b) k 
(c + n) k k\ 



s(n) := 



r(c + n) T(c + n-a-b) 
T(c + n- a) T(c + n-b) ' 



(2.3) 



The sum on the left-hand side of (2.2) still absolutely converges (since n > 0). Formula 
(2.2) can be rewritten as 



X>(n,fc) = l, 



(2.4) 



k=0 
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where 

F( k) = f{n ' k) = T(a + k) T(b + k)T(c + n- a) T(c + n- b) 
{U ' h s(n) F(l + k)T(c + n + k)r(a)r(b)r(c + n-a-b)' { } 

By (1.1) we have 

F(n, k) ~ fca+b-c-i-n as k _^ QO _ 

Thus the sum on the left-hand side of (2.4) converges absolutely (as we already knew). We 
want to prove (2.4) by the WZ method. Gosper's algorithm applied to F(n + 1, k) —F(n, k) 
or, equivalently, the WZ method applied to F(n, k) succeeds. In Maple V4, for instance, 
call 

read ~ hsum . mpl " : gosper(F(n+l ,k)-F(n,k) ,k) ; 

or 

read "hsum.mpl" : WZcertif icate(F(n,k) ,k,n) ; 
from Koepf's package hsum.mpl [4], or call 

read ekhad; ct (F(n,k) , 1 ,k,n,N) ; 
from Zeilberger's package EKHAD [8]. The resulting formula is: 

F(n + 1, k) - F(n, k) = G(n, k + 1) - G(n, k), (2.6) 

where 

G(n, k) := F(n, k). 

v ' c+n-a-b v ' 

Note that 

ua-\-b—c—n 

G(n,0)=0, G(n,k) as k -> oo. (2.7) 

v ' ' c+n-a-b v ' 

It follows from (2.6) and (2.7) that 

K K 

F ( n + 1 , k )~Yl F ( n > fc ) = G ( n > K + 1 )~ G ( n ' °) = K + l)^Q as K -> oo. 

k=0 k=0 

Hence X^fcLo F ( n ' ^) * s independent of n. Thus 

oo oo oo oo 

F(n, k) = lim ^ F(n, k) = £ ( lim F(n, k)) = ^ 5 k , = 1, (2.8) 

k=0 k=0 k=0 k=0 

where the second equality is justified by dominated convergence. Indeed, for some C > 
we have l/\s(n)\ < C if n G {0,1,2,...} (use (2.3) and (1.1)). Thus, by (2.5) and (2.3) we 
obtain 

|F(n,A0|<C|/(n,A0|< { " )k{h)k 



(2 9) 

(Rec + n ) k k\ 

for n G {no, no + 1,. . .}, where no is such that Rec + n > 0. Again by (1.1), the right- 
hand side of (2.9) summed over k = 0, 1, 2, . . . yields a convergent sum, thus justifying the 
application of the dominated convergence theorem in (2.8). 
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3. Non-terminating analogue of Saalschiitz' summation formula 

Saalschiitz' summation formula (which actually goes back to Pfaff, 1797) for a terminating 
Saalschiitzian 3 F 2 series with argument 1 reads as follows: 



3-^2 



— m, b, c 
e, — m + b + c — e + 1 

__ (e - b) m (e - c) m 
(e) m (e - b - c) m 



1 



E 



(-m)fc (b)k (c)k 



(e)k {—m + b + c - e + l) k k\ 



k=0 ' ' " (3.1) 
(m E {0, 1,2,.. .}, e, e - 6 - c £ {0, -1, -2, ...}). 



See for instance [1, 4.4(3) with proof in §2.1.5]. Now replace — m by a on the left-hand 
side of (3.1), where a is arbitrarily complex: 



3^2 



a, 6, c 
e,a + 6 + c — e + 1 



(a)fc (6)fc (c)fc 



fc=0 



(e) k (a + b + c - e + l) k k\ 



(3.2) 



and try to generalize identity (3.1) for this case. Such generalizations, with one additional 
3^(1) term, are known in literature (see for instance ...), but we want to find a formula of 
this type from scratch, just starting from (3.2) and working with the Zeilberger algorithm. 
In (3.2) replace c by c + n, where n G {0, 1,2,.. .}. Then 



3^2 



a, 6, c + n 
e,a + 6 + c — e + n + 1 



E /fa*), 



(3.3) 



k=0 



where 



so 



f(n, k) 



(a)k (b)k (c + n) k 



(e) k (a + b + c-e + l + n) k k\ 



r e T (a + b + c-e + l + n , _ 2 
7 r(a)T(6)r(c + n) 

by (1.1). Thus we have absolute convergence in the sum in (3.3). Zeilberger's algorithm 
applied to /(n, k) succeeds. In Maple V4, for instance, call 

read ekhad: ct (f (n,k) , 1 ,k,n,N) ; 

from Zeilberger's package EKHAD [8]. The resulting recurrence is: 

(c - e + a + n + 1) (b + c - e + n + 1) f(n + 1, k) 

-(c-e + n + l)(b + c-e + a + n + l) f(n, k) = g(n, k + 1) - g(n, k), (3.4) 



where 



, (b + c-e + a + n + l)(e + k-l)k r . 

g{n, k) := ■ f(n, k). 

c + n 
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We can rewrite (3.4) in the form 

F(n + l,k)- F(n, k) = G(n, k + 1) - G(n, k), (3.5) 

where 

„, , % T(a + c-e + l + n)T(b + c-e + l + n) 

F(n,k) :=— -. — — - — ^ - f(n,k), 

v ' 1 F(c-e + l + n)r(a + b + c-e + l + n) JK ' h 

Observe that 

r(e)r(a + c-e + l + n)r(6 + c-e + l + n) ,_ 2 
F(n, k) ~ — — — — 77^ — / ; k as k -> oo, 



and 

rv m n r r( 1A r(e) r(g + c - e + 1 + n) T(b + c - e + 1 + n) 

G (n, = 0, lim G (n, k = -- —77—— — — r . (3.6 

v ' 1 ' fc^oo v ' 1 r(a)r(6)r(c+l + n)r(c-e + 2 + n) v ; 

Put 

00 

S(n) :=J2F(n,k) 



fc=0 

^T(a + c-e + l + n)T(b + c-e + l + n) 
~~ r(c - e + 1 + n) T(a + b + c - e + 1 + n) 3 2 

Then by (3.5) and (3.6) we obtain 



a, 6, c + n 
e, a + 6 + c — e + l + n' 



• (3-7) 



S(n + 1) - S(n) = lim G(n, k + 1) - G(n, 0) 

fc^oo 

= r(e) r(a + c - e + 1 + n) Tip + c - e + 1 + n) 

r(a)r(6)r(c + i + n)r(c-e + 2 + ?i) ' ^ ' ' 

Hence 

q( , _ q( . T(e)T(a + c-e + l)T(b + c-e + l) ^ (a + c - e + l) fc (6 + c - e + l) k 

[n) u r(o) r(6) r( c + 1) r( c - e + 2) ^ ( c + i) k ( c - e + 2), 

(3.9) 

Note that 

(g + c-e + l) fc (6 + c-e + l) fc T(c + 1) T(c - e + 2) , a + 6 - e -i nc , _ 

(c + l) fc (c - e + 2) fc T(a + c - e + 1) T(6 + c - e + 1) 
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Therefore assume that Re (e — a — b) > 0. Then, for n — > oo, the right-hand side of (3.9) 
tends to 



T(a + c - e + l)T(b + c - e + 1) 
r(c - e + 1) T(o + b + c - e + 1) 



3-^2 



a, 6, c 



;i 



r(e) r(a + c - e + 1) T(b + c - e + 1) 

r(a)r(6)r(c + i)r( c -e + 2) 3 2 



e, a + 6 + c — e + 1 

a + c — e + l,b + c — e + 1,1 



c+l,c-e + 2 



;i 



Here we substituted the right-hand side of (3.7) with n = for 5(0). Next let n — > oo on 
the left-hand side of (3.9). From (3.7) we obtain 



lim S(n) = iF\ 



a, b 



r(e)r(e-a-ft) 
r(e-a)r(e-6) 



(Re(e-a-ft) > 0). 



(3.10) 



The second identity is Gauss' summation formula (2.1). The first identity follows by 
applying (1.1) to the limit of the quotient of Gamma functions in (3.7) combined with a 
formal limit (to be justified in the Lemma below) within the 3^2(1) m (3-7). Thus the 
limit case of (3.9) for n — > 00 is the identity we looked for: 



T(a + c - e + 1) r(6 + c - e + 1) 
T(c - e + 1) T(a + b + c - e + 1) 



3-^2 



r(e) T(a + c - e + 1) T(6 + c - e + 1) 

+ r(a)r(6)r(c + i)r( c -e + 2) 3 2 

_ r(e)r(e-a-6) 
~ r(e - a) r(e - 6) ' 



a, 6, c 
e, a + b + c — e + 1 ' 

a + c — e + 1, 6 + c — e + 1,1 



c+l,c-e + 2 



;i 

(3.11) 



valid for Re (e — a — b) > 0. 

Let us now give the promised Lemma which will justify the first identity in (3.10). 



Lemma 3.1 Let Re (e — a — b) > 0. Then 



lim 3 F 2 



a,b,c + n 
e,a + b + c — e + 1 + n 



;i 



— 2F1 



a,b 
e 



Proof The limit formally holds because 

(a)k (b)k (c + n) k 



lim 



(a)k (b)t 



00 (e)k (a + b + c - e + 1 + n) k k\ (e) k k\ 

We will justify this limit by dominated convergence. Take e £ (0, 1) such that also e < 
Re (e — a — b). It follows from (1.1) that 
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(a)k (b)k (c + n) k 



const. 



(e)k (a + b + c - e + 1 + n) k k\ 

T(c + a + b - e + 1 + n) 



x 



T(c + n) 



T(a + k)T(b + k) 
r(e + k) r(i + /c) 

r(c + n + fc) 



r(c + a + 6- e + l + n + A;) 

< Const, fcRe(a+6-e)-l n Re(a+&-e) + l ^ + fc jRe(e-a-6)-l 

(n" 1 + fc -l)Re(e-a-6)- £ 



= const. 



l-e 



n + k 



k 



-l-e 



( Const. fc-Re(e-a-6)-l < con st. /c 

I const. k~ e ~ l 



- - 1 if fc < n, 
if /c > n. 



Then the dominated convergence follows because YlkLi ^ £ 1 



< oo. 



□ 



4. Other non-terminating identities related to Saalschiitz' summation formula 

Formula (3.11), which I derived in the previous section by a mixture of computer algebra 
and asymptotic techniques, is possibly not in the literature, but it can be derived from 
some formulas of similar nature which are in the literature. 

The best known non-terminating generalization of Saalschiitz' formula (3.1) is: 

a, 6, c 



3 r 2 



e, a + b + c — e + 1 ' 

r(e - 1) T(a - e + 1) T(6 - e + 1) T(c - e + 1) T(a + 6 + c - e + 1) 



+ 



r(l - e) T(a) r(6) r(c) r(a + b + c - 2e + 2) 

a — e + 1, b — e + 1, c — e + 1 
2-e,a + 6 + c-2e + 2 ' 

T(a - e + 1) r(6 - e + 1) T(c - e + 1) T(a + b + c - e + 1) 



(4.1) 



T(l - e) r(6 + c - e + 1) T(a + c - e + 1) T(a + b - e + 1) 

see [2, (11.24)]. Since both 3^(1) series are Saalschiitzian, there are no conditions on the 
parameters necessary for convergence. 

The second formula we will use involves two non-terminating 3^2(1) series with one 
upper parameter equal to 1, like the second 3-F 2 (l) in formula (3.11). The formula reads 
as follows. 



3^2 



a, 6, 1 
cL e 



;1 



+ 



l-e 
e-b-l 



3^2 



d — a, b, 1 
d, b - e + 2 : 



1 



T(d) r(e) r(a - e + 1) T(6 - e + 1) + e - a - & - 1) 

r(o)r(6)r(d-a)r(d-6) 

(Re (d + e - a - 6 - 1) > 0, Re (a - e + 1) > 0). 



(4.2) 
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It can be derived by specialization of [2, (3.3.1)] combined with Gauss' summation formula 
(2-1). 

The third formula transforms a nonterminating Saalschiitzian 3^2(1) into a nonter- 
minating 3-F 2 (l) with upper parameter 1: 



3-^2 



a, 6, c 
e, a + b + c — e + 1 



;i 



r(e)r(a + 6 + c-e + l) 

r(a)r(6 + i)r( c + i) 



X 3 ^2 



e — a,b + c — e+1,1 
b+ l,c+ 1 



;i 



(Rea>0). (4.3) 



It follows by specialization of [2, (3.1.2)]. 

Now replace in (4.3) a, 6, c, e by 6 — e + 1, a — e + 1, c — e + 1, 2 — e. Next substitute in 
the right-hand side of (4.3) (with parameters replaced as above) formula (4.2) with a, 6, d, e 
replaced by 1 — 6, a + c — e + 1, a — e + 2,c — e + 2. This yields 



3-^2 



6 — e + l,a — e+l,c — e + 1 



2-e, a + 6 + c-2e + 2 



X3.F2 



r(2-e)r(a + 6 + c-2e + 2) 
" cT(6 - e + 1) T(a - e + 1) T(c - e + 2) 



1 + b + c — e, 1 + a + c — e,l 



2 + c-e, c+1 

r(e - a - b) r(c) r(2 - e) T(a + 6 + c - 2e + 2) 
"r(l - 6) T(l + 6 + c - e) T(l - a) T(a + c - e + 1) 

(Re (1 + b - e) > 0, Re (e - a - 6) > 0) 



(4.4) 



Now formula (3.11) follows from (4.1) by substituting (4.4) for the second 3-F 2 (l) in (4.1), 
and by using that 



Y{z)Y{l-z) 



n 



sm(nz) 



So we have seen that formulas (4.1), (4.2) and (4.3) imply formula (3.11) via formula 
(4.4). Conversely, formulas (3.11), (4.2) and (4.3) imply formula (4.1) via formula (4.4). 
In the next section we will derive formulas (4.2) and (4.3) by methods of computer algebra 
and asymptotics. 



5. Further proofs of identities of non-terminating series by computer algebra 
and asymptotics 

Let us try to prove formula (4.2) in a similar way as formula (3.11). Put 

f(n,k):= ff» ( V7t . 

(d + n) k (e)fc 
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By Zeilberger's algorithm we find that 

(a - n - d) (n + b) f(n + 1, k) - (n + b - e + 1) (n + d) f(n, k) = g(n, k + 1) - g(n, k), 

where g(n, k) := —{n + d) (e + k — 1) /(n, k). 



Hence 
where 



F(n,k) := 
G(n,k) := 



F{n + 1, k) - F(n, k) = G(n, k + 1) - G(n, k), 
(b) n (d-a) n (a) k (b + n) k 



Assume 



Then 



(d) n (b-e + l) n (d + n) k (e) k ' 

(e + k-l)(a) k (b + n) k (b) n (d-a) n 
(e - n - b - 1) (d + n) k (e) k (d) n (b - e + l) r 

Re(d+e-a-b- 1) > 0, Re (a - e + 1) > 0. 



S(n) :=J2F(n,k) 



k=0 



(b) n (d-a) 
(d) n (b-e+l) n 



3-^2 



a, b + n, 1 
<i + n, e 



has absolutely convergent series and 



S(n + 1) - 5(n) = lim G(n, k + 1)- G{n, 0) = -G(n, 0) 

fc^oo 

(e - 1) (d - a) n (b) n 



(e-6-l)(d) n (6-e + 2), 



Hence 



Thus 



S(0 ) = S(oo) + __ g (<i)n(6 _ e + 2)ii , 

where S'(oo) := lim S(n). 

n^oo 

e — 1 „ |~ <i — a, 6, 1 



3-^2 



a, 6, 1 
<i, e 



;i 



e-b-l 



3^2 



d, 6 - e + 2 



;i 



+ 5(oo). 



So formula (4.2) will be proved if we can show that 

a, b + n, 1 



lim (6) n (d-q)n ^ 

n^oo (d) n (b - e + l) n 



d + n,e 



;i 



T(d) r(e) r(a - e + 1) T(6 - e + 1) T(d + e - a - b - 1) 

r(a)r(6)r(d-o)r(d-6) 



(5.1) 
(5.2) 



(5.3) 



(5.4) 
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This limit result is not evident. A formal limit would yield 

o, 1 



( lim n 6 ""- 1 ) 2 Fi 



;i 



= 0-oo since Re (a — e + 1) > 0. 



We will give a proof of (5.4) at the end of this section. 

Let us first turn to formula (4.3), and rewrite it in the following equivalent form. 



(b) n (d-a) n 
(d) n (6-e+l), 



3-^2 



a, b + n, 1 
d + n,e 



;i 



r(d) r(e) r(6 - e + 1) T(d + e - a - b - 1) 



T(b + n)T(d - a + n) 



r(6) T(d - a) T(d + e - b - 1) r(6 - e + 1 + n) T(d + e - a - 1 + n) 

d + n — l,d + e — a — b — 1, e — 1 



d + e — a + n — l,d + e — b — 1 



(Re(d+e-a-6-l) >0). (5.5) 



Observe that we can take limits on the right-hand side of (5.5) for n — > oo (see (3.10)), 
and thus arrive at (5.4). However, I only want to use (5.5) if I can derive it by methods of 
computer algebra and asymptotics. Let us try. 
Write (5.5) more compactly as 



S(n) =CS(n), 

where S(n) is given by (5.2), 

T(d) r(e) r(6 - e + 1) T(d + e - a - b - 1) 



(5.6) 



C 



T{b)Y{d - a)T(d+ e - b - 1) 



and 

S(n) 



T(b + n)T(d - a + n) 



T(b-e + l + n)T(d + e~a-l + n) 



3-^2 



d + n — l,d + e — a — b — 1, e — 1 
d + e — a + n — l,d + e — b — 1 



;1 



Note that <S(n) is just (3.7) with a, 6, c, e replaced by e — 1, d+e — a — b— 1, <i— 1, d+e — b — 1. 
Thus (3.8) yields that 



S(n+1)-S(n) 



1 (e - 1) (d - a) n (6), 



C (6-e+l)(d) n (6-e + 2) n ' 
Comparison with (5.3) yields that S(n) — C S(n) is independent of n. Hence 

S(n) -CS(n) = S(oo) -CS(oo). 

Thus (5.6) holds iff (5.4) holds. So we have not made real progress. Formula (5.4) will follow 
from (5.6), but for the proof of (5.6) we need (5.4). Let us therefore give an independent 
proof of (5.4). 
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Proof of (5.4) Observe that 

a, b + n, 1 



(b) n (d-a) n 

3-^2 



;i 

a + n, e 



(d)„(6-e + l) n 

_ T(d) r(e) r(& - e + 1) y> T(a + fc) r(6 + n + fc) r(rf - a + n) 

7 " r(a)T(6)r(d-a) ^ r(e + fc) r(d + n + k) T(b-e + l + n) ' ^ ' ' 

We will first work formally, and consider the sum obtained by replacing in the above 
infinite sum the gamma quotients by their asymptotic estimates (for big n and fc) as given 
by (1.1). This yields 

y> T(q + fc) T(b + n + fc) T(d-a + n) 

^ r(e + fc) T(d + n + k) T(b - e + 1 + n) ^ ' ' 

X) 

J2(k + l) a_e (n + k + 2) b ~ d (n + if+e-a-b-i ^ 

pOC 

■ / (1 + t) b - d t a ~ e dt asn^oo (5.11) 

Jt=0 

T(a - e + 1) r(d + e - a - b - 1) 

~ r(d-b) ' 

Here we used that the integral in (5.11), which converges because of assumptions (5.1), is 
approximated by its Riemann sum (5.10) asn^ oo. Thus formally the right-hand side of 
(5.7) indeed equals the right-hand side of (5.4). 

We will justify this formal proof by Lebesgue's dominated convergence theorem. 
Rewrite (5.8) as 

pOO 

ln {t)h n {t)dt, (5.12) 



/ 

Jo 



where, for ^ < t < 

K{t) := ( 1 + 



k + l\ b ~ d /fc+1 



n + 1 J \n + 1 
and 

T(a + k) T(b + n + k) T(d-a + n) 

7n ^ "~ r(e + fc) (fc + l) a " e T(d + n + k) (n + k + 2) h ~ d T(b - e + 1 + n) (n + i)d+e-a-b-i 
Then 

hn(t)\ < const, and |/i n (t)| < const. (1 + t) Re(b " d) t Re(a " e) . 

Furthermore, 

lim ln (t) = 1 and lim h n (t) = (1 + t) b ~ d t a ~ e . 

n^oo n— >oo 

It follows by dominated convergence that the integral (5.12) tends to the integral (5.11) 
as n — > oo. EH 
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